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Abstract 

We study here the structure of singularity forming in gravitational collapse of spher¬ 
ically symmetric inhomogeneous dust. Such a collapse is described by the Tolman- 
Bondi-Lemaitre metric, which is a two-parameter family of solutions to Einstein equa¬ 
tions, characterized by two free functions of the radial coordinate, namely the ‘mass 
function’ F(r) and the ‘energy function’ f{r). The main new result here relates, in 
a general way, the formation of black holes and naked shell-focusing singularities re¬ 
sulting as the final fate of such a collapse to the generic form of regular initial data. 
Such a data is characterized in terms of the density and velocity profiles of the matter, 
specihed on an initial time slice from which the collapse commences. Several issues 
regarding the strength and stability of these singularities, when they are naked, are 
examined with the help of the analysis developed here. In particular, it is seen that 
strong curvature naked singularities can develop from a generic form of initial data in 
terms of the initial density profiles for the collapsing configuration. We also establish 
here that similar results hold for black hole formation as well. We also discuss here the 
physical constraints on the initial data for avoiding shell-crossing singularities; and also 
the shell-focusing naked singularities, so that the collapse will necessarily end as a black 
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hole, preserving the cosmic censorship. These results generalize several earlier works 
on inhomogeneous dust collapse as special cases, and provide a clearer insight into the 
phenomena of black hole and naked singularity formation in gravitational collapse. 
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I. Introduction 


The classic paper by Oppenheimer and Snyder analyzes the gravitational collapse of 
a spherically symmetric massive object to conclnde that the collapse of a homogeneous dust 
cloud would result into a black hole in the spacetime. Such a black hole is characterized by 
the presence of an event horizon, and the spacetime singularity at the center, which is covered 
by the horizon. This scenario provides the basic motivation for the black hole physics, and 
the cosmic censorship conjecture 0, which states that even when the assumptions contained 
in the above case are relaxed, in the form of either perturbations in the symmetry or form 
of matter etc., the outcome would still be a black hole in generic situations. As no proof, or 
a rigorous mathematical formulation of the censorship hypothesis has been available so far, 
the course of action that has emerged in past decade or so has been to carry out a detailed 
investigation of different gravitational collapse scenarios in general relativity to understand 
the final fate of collapse ( see e.g. and references therein). 

Within such a context, an important situation that immediately offers itself for analysis is 
the introduction of inhomogeneities in the matter distribution. It is clear that the assumption 
of homogeneity is only an idealization, and realistic density prohles, for massive objects such 
as stars will have inhomogeneous density distribution, peaked typically at the center of 
the object. One would like to examine the final fate of gravitational collapse of such an 
inhomogeneous dust cloud, and examine in what ways there are departures in conclusions 
as opposed to the homogeneous situation. The Einstein equations can be fully solved for 
the stress-energy tensor of the form of inhomogeneous dust, and such a collapse is described 
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by the Tolman-Bondi-Lemaitre (TBL) metric. This is a two-parameter family of solutions, 
characterized by two free functions of the radial coordinate which are the ‘mass function’ 
F{r) and the ‘energy function’ /(r) of the cloud. While the hrst quantity describes here 
the initial density distribution for the collapsing cloud, the second function characterizes the 
initial velocity prohle of the collapsing shells. 

Our main purpose here is to relate, in a general way, the formation of black holes and 
naked shell-focusing singularities resulting as the final fate of such a collapse to the generic 
form of regular initial data. Such a data is characterized in terms of the density and velocity 
prohles of the matter, as given by the two functions above, specihed on an initial time 
slice from which the collapse commences. An important issue, when naked singularities do 
form in gravitational collapse, is that of their curvature strength. We discuss here several 
issues regarding the strength and stability of such singularities, with the help of the analysis 
developed here. In particular, it is seen that strong curvature naked singularities can develop 
from a generic form of initial data in terms of the initial density prohles for the collapsing 
conhguration. We also establish here that similar results hold for black hole formation 
as well. Thus, we show that both the black hole as well as naked singularity formation 
are related to the nature of the regular initial data dehned at the onset of gravitational 
collapse. We also discuss here the physical constraints on initial data for avoiding shell¬ 
crossing singularities; and also the shell-focusing naked singularities, so that the collapse will 
necessarily end as a black hole in conformity with the cosmic censorship conjecture. These 
results generalize several earlier works on inhomogeneous dust collapse as special cases. 
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and provide a clearer insight into the phenomena of the black hole and naked singnlarity 


formation in gravitational collapse. For example, while some of the earlier works ( see e.g. 
05 0 0’ 0) discnss the formation of naked singnlarities in dnst collapse, either nnder varions 
special assnmptions or nnder general conditions, they do not address the important physical 
problem of relating this naked singnlarity formation to the natnre of regnlar initial data 
at the onset of the gravitational collapse. It is crncial, from the physical point of view, 
to characterize and identify the properties of initial data which distingnish between the 
formation of black holes and the occurrence of naked singularity. Such an analysis could also 
pave the way for any possible formulation and proof for the cosmic censorship conjecture. 
A beginning in such a direction was made in We present here a more general and 

complete analysis, obtaining several new results in connection with the naked singularity 
and black hole formation in gravitational collapse, in the specihc context of the genericity of 
the initial data. 

The main issue of concern here to us is, given a generic density distribution at the onset of 
collapse what all is possible as the hnal outcome. It is seen that given such a generic density 
prohle, one can always have either a black hole or a strong curvature naked singularity, 
depending on the nature of the initial data. Certain points related to genericity are also 
discussed. 

Section II below provides a brief review of the TBL models, giving the main set of 
equations used. Section III specihes the generic initial data to be used for collapse, and 
discusses the formation of singularities in collapse situations. Shell-crossing singularities are 
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discussed in Section IV. The formation of black holes and naked singularities is discussed in 


Section V, and the last section gives certain concluding remarks on our results. 


II. Tolman-Bondi-Lemaitre Dust Models 


The TBL spacetimes are spherically symmetric manifolds (M,g), with metric of the form, 


ds'^ = —dt^ + e^'^dr"^ + R^dVP', 


( 1 ) 


and energy-momentum tensor of the form of a perfect fluid with equation of state p=0, 
given by 


= €6161 


( 2 ) 


Here e, u and R are functions of r and f, and is the metric on the 2-sphere. The 
Einstein equations become 



(3) 

(1 + /(r)) = 

(4) 

. X F'(r) 

(5) 
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where the dot and prime signify partial derivatives with respect to t and r respectively. 
It is seen that the energy density blows np either at the “shell-focnsing singnlarity” R = 0, 
or when R' = 0 which corresponds to a “shell-crossing singnlarity” in the spacetime. The 
fnnctions F{r) and /(r) are free fnnctions of integration which completely specify the initial 
data in the model. Here F{r) can be interpreted as the weighted mass given by 


F{r) = [ {1 + f)e{r,t)dvt = 4:71 [ pr'^dr 
JBr-A Jo 


( 6 ) 


where p = e(0,r), and is a ball of coordinate radins r centered on r=0 in the hyper- 
snrface = const.). The time slicing (t=const.) has been so chosen snch that r=const. 
labels a matter shell with t measnring the proper time elapsed along its geodesic path, that 
is we write the metric in comoving coordinates. Eqnation (3) can be thonght of as the 
general relativistic generalization of the Newtonian energy eqnation which leads to the 
interpretation of F{r) as the effective mass contained in the sphere of radius r; and /(r) 


as effective total energy per unit mass of a fluid element labeled by the comoving radial 
coordinate r. The model is said to be bound, unbound, or marginally bound if /(r) is less 
than zero, greater than zero, or equal to zero respectively. 


The integrated form of equation (3) is given by 


t - tJr) = 


R^/^G{-fR/F) 

7f 


(7) 


where G{y) is a real, positive and smooth function which is bounded, monotonically 
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increasing and strictly convex, and is given by 


r arcsin^ _ yi-y 

y3/2 y ) 


1 > 1 / > 0 , 


G{y) = I 


y = 0, 


-arcsinhy/^^ _ (1 > 7/ > -OO 

I (-27)3/2 y , U > 7/ ^ CX) 


( 8 ) 


and ts{r) is a constant of integration which can be hxed by the choice of scaling on the 


initial surface {t = 0). Using this scaling freedom, if we choose so that i?(0,r) = r then. 


isir-) 


r^l‘^G{-fr/F) 

7f 


(9) 


The time t = ts{r) corresponds to the value i? = 0 where the area of the matter shell at 
a constant value of the coordinate r vanishes, which corresponds to the physical spacetime 
singularity. Thus the range of coordinates is given by 


0 < r < Tc, —CX) < t < ts{r), 


( 10 ) 


where r = Vc denotes the boundary of the dust cloud where the solution is matched to the 
exterior Schwarzschild solution. 

The argument of the function G{—fR/F) changes it’s sign depending on the sign of f{r), 
since both R and F{r) take only non-negative values. If /(r) > 0 for all r > 0, then from 
equation (3) we get i? 7 ^ 0 for all t. From the integrated form of equation (3) we can show 
that if for some r > 0 we have /(r) > 0 then this particular shell is strictly unbounded, and 



from a physical point of view, the collapse can be interpreted as not being gravitational, but 
due to specihc initial conditions (inhnitely far in the past) [^. Therefore, we shall discuss 
here only the more physical case with /(r) < 0 (however, it will be possible to generalize the 
conclusions here to other cases as well by similar methods as used here). 

It is seen that if the condition 

R' >0 (11) 

is satished then the dust collapse always avoids shell-crossing singularities. The restrictions 
implied by this condition on the initial data will be discussed in sections IV and V. The 
weak energy condition, i.e. Rj = e{rR)V^V^ > 0, for all nonspacelike vectors, V* is assumed 
everywhere in spacetime. This imposes a condition on the initial data that the energy density 
e(r, t) is non-negative everywhere. 

III. Generic Initial Data and Formation of Singularities 


We call a spacetime to be singular [^, p!2 
geodesic. 


if it contains an incomplete nonspacelike 


7 : [0, x) ^ At 

(where Ai is the spacetime manifold) such that there is no extension 


9 






Q:M^ M! 


for which ^07 is extendible. The spacetime should be inextendible, that is, we should not 
be able to extend the above incomplete curve by embedding the spacetime into a larger man¬ 
ifold. This essentially means that we do not admit singularities created by removing pieces 
from the spacetime manifold. Important theorems showing the existence of singularities in a 
spacetime under “reasonable” physical conditions, which might form either in gravitational 
collapse or cosmology, were proved by Penrose, Hawking, and Geroch jl^ , |T^. While these 
theorems prove the existence of singularities under a fairly broad spacetime framework, they 
provide no information on the nature and structure of the singularities they predict, either 
in collapse scenarios or cosmology. In particular, they do not imply whether the singularities 
resulting as the final fate of gravitational collapse will be necessarily covered, or otherwise, 
by the event horizons of gravity. Our purpose here is to analyse the nature and structure 
of curvature singularities occurring in the TBL collapse models. A singularity will be called 
naked if there are future directed nonspacelike curves in the spacetime which terminate at 
the singularity in the past, otherwise it will be hidden behind a black hole. Considering 
equation (7), the apparent horizon within the collapsing dust cloud lies at i? = F{r), i.e. 

tahir)=tsir)-FG{-f) (12) 

It can be easily seen from the above equation that all the points on the singularity 
curve ts{r), other than the central point (r = 0) are covered by the apparent horizon. This is 
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because, since both the functions F{r) and G{r) are strictly positive for r > 0, with F{r) = 0 
at r = 0, therefore for all r > 0 ts(r) > tah{r) and fs(0) = fa/i(0). Thus, all the other points 
on the singularity curve, except the point r = 0, are covered by the apparent horizon. It is 
the central singularity point r = 0 whose nature, in terms of being naked or covered, depends 
on the data specihed at the initial surface. 

In the TBL models the initial data, to be specihed on the initial hypersurface Ej, consists 
of two free functions F{r) and /(r) respectively and such a specihcation uniquely determines 
the solution to held equations inside which is the future Cauchy development of the 

hypersurface E*. To study the formation of singularity for this initial data, we consider the 
initial functions in a “generic” (most general expandable) form 

OO OO 

Fir) = T. Fnr''+^ S{r) = /„r" (13) 

n=0 n=2 

The choice of n here is made to avoid any singular behaviour of functions, and to ensure 
the hniteness of density, on the initial surface. The functions chosen here are assumed to be 
only, and not necessarily analytic (C"^); since the only requirement of the theory is to have 
functions rather than restricting to the stronger diherentiability condition of analyticity. 
This way, we are able to deal with a broader class of spacetimes. Additional differentiability 
requirements higher than C‘^ will make the situation functionally less generic. We now recall 
briefly some dehnitions [0 which are necessary for further analysis. Especially, the partial 
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derivative for the area function R is written as, 


R' = r“-i[(r/ - f3)X + [e-{ri- lP)X^/^G{-PX)] x [P + ^]V2] 
= r“-i/7(X,r) 


where we have used the following 


(14) 


F' f f 

u = r°‘, X = {R/u), rjir-) = r—, /5(r) = r^, p(r) = r^, P(r) = pr^~^, 

F f F 

A = £ 1 + P-ri (r,-lP)G(-p) 

^a-1 _|_ p^l/2^3(a-l)/2 ^3(a-l)/2 

The function P{r) is dehned to be zero when /(r) is zero (the marginally bound case). 
The parameter a (which satishes a > 1) is introduced here for examining the structure of 
the central singularity at r = 0. In terms of these new variables, we can write for energy 
density as 


rjA 

pm 


( 16 ) 


The weak energy condition implies H{X,r) > 0, and pA > 0. The actual value of a is 
uniquely determined by the initial data prescribed, and is the key factor in deciding whether 
or not we have a strong curvature (covered or naked) singularity. The expression for energy 
density on the initial surface (which is the scaling surface Sj, on which R = r) is e = F'/r^, 
and the weak energy condition implies that F' > 0. The requirement that initial surface 
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should not contain any trapped surfaces {F{r) /R{r,t) > 1) also gives a restriction on the 
choice of initial data. 

The equation for radial null geodesics in the spacetime (1) is given by 


M _ R' 

dr (1 + f{r)yR 


(17) 


This can also be written in the form, 


^ = n _ + 

du a/1 + / a 


U{X,r) 


( 18 ) 


Dehning 


Xo = lim — = lim (19) 

R^o u R^o du 

w —>0 hi —>0 

it can be shown that if the equation 

a//o + Aq / X 

V{X) = U{X, 0) - X = (1 - 

admits a real positive root, then the central singularity at r = 0, i? = 0 is naked (at least 
locally). The global visibility of such a singularity will depend on the overall behaviour of 
the concerned functions within the dust cloud in the range 0 < r < r,,. In the case otherwise, 
a black hole will be formed as the end product of collapse. The parameter a is uniquely 


H{X,0) 


X = 0 


a 


( 20 ) 
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fixed by demanding that 0(X, r) goes to a nonzero finite value in the limit r —0. Consider 
the expression for 0 


l + f3-ri (r/- |/3)G'(-p) _ 'h(r) 

~ ^a-l ~ _|_ p^1/2^3(q- 1)/2 ^3(a-l)/2 ~ ^3(a-l)/2 


( 21 ) 


where 


\l/(r) = 'ho + + ^2r‘^ + + ^^4/ H- 


( 22 ) 


The value of a is decided by the hrst nonvanishing term in the expansion of 'l/(r). Each 
term in the expansion of 'h(r) is completely specihed in terms of the initial data; so given 
such data at the initial surface in terms of the density and velocity distributions, as specihed 
by the functions F and / above, we can tell using the criterion above whether the hnal fate 
of collapse is going to be a black hole or a naked singularity. 


The nonspacelike geodesics of the spacetime (1) are given by. 




dt 

dk 


V 

R 


^ _ ^/TTjyV^ -P + BR^ 
dk RR' 


{K^f + sin^e{K^f = f/R"^ 


(23) 


(24) 


(25) 


where the = dx^'/dk denote the tangent to the geodesics. Here k is affine parameter 
along geodesics, I is the impact parameter (/ = 0 ^ radial trajectories), and the values 
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B = 0,-1 correspond to null and timelike curves respectively. The function V satisfies the 
differential equation 


dV 

dk 


+ (p2 


l^ + BR^)[ 


R' 

RR 


R^ R^ 


(26) 


For the clarity of discussion, we consider only radial null geodesics here, in which case 
the above equation reduces to 


1 dV R R VTTJ 
lp‘ IP 


(27) 


We shall call the singularities to be strong (see e.g. 
tories the following is satished. 


1 ^) if along nonspacelike trajec- 


lim k‘^'i/j{r) = lim k'^RijV^V^ ^ 0 (28) 

fc —>0 fc —>0 

Such a rate of divergence indicates a very powerful curvature growth in the limit of 
approach to the singularity (which is the same as in the case of the big bang singularity of 
cosmology), and ensures that all the volume forms, as dehned by the Jacobi helds along the 
nonspacelike geodesics, are crushed to zero size in the limit of approach to the singularity. 
Physically, this can be interpreted as indicating that all the objects falling into the singularity 
are crushed to zero size. The signihcance of this would be that there could not possibly be 
any extension of the spacetime through such a singularity, as opposed to a gravitationally 
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weak singularity such as a shell-cross through which the spacetime could possibly be extended 
and continued further. 

In case of TBL models, we have 


F'( 

\imk‘^^jJ(r) = lim-- —- 


(29) 


For radial null geodesics, we can write using the expressions above for R, R', and using 
the L’Hospital rule. 


lim 

k^O 


y 4rioAoffo 

lim- — - 

k^o Xo^[Vl + / o ( 3 q ( — i]o) ~ -^o]^ 


(30) 


Where the quantity Nq is the limiting value of —rR', and Aq is given by 


'0, 

q; < 3, 


i"o, 

0 = 3, 

(31) 

. CX), 

q; > 3. 



Since Aq occurs in the numerator of expression for k'^'ipijr), we have strong curvature 
singularities only for a > 3. All the same, it follows from equation (18) that whenever a 
is greater than 3 we always have black hole, since we cannot have any outgoing geodesics 
meeting the singularity in the past in that case. 


IV. Shell-Crossing Singularities 
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Some counterexamples to cosmic censorship were proposed [119 using the so called shell 


crossing singularities by Yodzis et. ah, who showed the existence of such singularities and 
also that they are naked. However, these singularities are gravitationally weak (both in 
Tipler as well as the Krolak |]Y] sense); and there have been proposals for extending the 


spacetime through such singularities, in particular, by Papapetrou and Hamui and also 
others. These singularities are generally not considered as being any serious counterexample 
to cosmic censorship conjecture, as the spacetime may be extended in a distributional sense 
through such a mild singularity. However, it is of importance to see what are the restrictions 
imposed on the initial data and the spacetime by the avoidance of such singularities at any 
point other then at the center, before the central shell-focusing singularity occurs; and to 
examine how ‘physical’ or ‘unphysical’ such restrictions are. That is, we would like to know 
what constraints the avoidance these singularities places on the functional form of initial 
data. We will impose the criterion for avoidance of shell-cross singularities on the initial 
data for making it physically more reasonable. 


Shell-crossings are singularities where we observe crossing of dust shells, and the density 
diverges at these points. Consider the expression for energy density (5) in the TBL models. 
We have a shell-crossing singularity at i?' = 0 and at this point there is a divergence in 
density and also certain curvatures. But for the spherically symmetric dust collapse it is 
possible to hnd a coordinate system so that we have a regular extension of the metric 
through these singularities, but which need not be 00. Calculating the expression for 
R' using equation (7), 
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[^-{^f/‘^G{sR/r) + qs{^fl‘^G\sR/r)]E! 


R n/2r3G{sR/r) m 

+g \ +p{-)G\sR r)] 

r 2 r 


r5s'[G^(s) — (—)^'^^G^(si^/r)] 
R[G{s) - {^fl^G(sRlr)] (32) 


where we have used notation, 


G^{x) = dG{x)/dx, s{r) = q{r) = (33) 

r [r) 


Clearly, if G{x) is strictly positive and convex, and R/r < 1 for f > 0, we cannot have 
shell crossings p if 


s' > 0 , g' < 0 


(34) 


The restrictions imposed by the above equations are physically reasonable in the following 
sense. The condition g' < 0 implies that the density should be constant or decreasing away 
from the center, which gives a restriction on the nonvanishing terms in the expansion for 
density. Now analysing the condition s' > 0, the kinetic energy term T and the potential 
energy term V (0, r) on the initial hypersurface are given by 


T = R\ 1/(0, r) 


F(r) 


(35) 
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Then from equation (3) we can write, 


T = -V{r) + /(r) 


(36) 


Using the above two expressions we can write, 


s(r) 


fjr) 

U(0,r) 


1 

1 - T/f{r) 


(37) 


The function s{r) takes values between 0 and 1. Since s'{r) > 0 (it is zero for homogeneous 
and marginally bound case), it is clear from the above equation that as we move away from 
the center to the boundary of the star, contribution of the kinetic energy term to total energy 
decreases. Physically this implies that we do not give additional velocity to the outer layers 
of the matter in the cloud, in comparison to the inner layers, to avoid shell crossings. 


V. Naked Singularities and Black Holes 


The singularity theorems prove only the existence, and do not give any information on 
the nature and behaviour of the singularities they predict. This leaves an important gap 
between the theoretical existence and physical presence of the singularities of cosmology 
and gravitational collapse. As mentioned earlier, one of the main differences lie in the issue 
related to genericity. As we discussed, in dust models the most general solution to Einstein 
equations contains two free functions, which are functions of the comoving radial coordinate 
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r. If the final outcome of gravitational collapse is specified in terms of both the functions 
(in a general form), the solution can be called as functionally generic (see e.g. also), at 
least within the context of the dust models we are discussing. 

We have discussed in Section III the conditions on the initial data for the formation of a 
naked singularity or a black hole. As was pointed out there, the behaviour of the hrst point 
r = 0 of the singularity curve depends on the initial data, and all other points on the curve 
are covered by the apparent horizon independently of the initial conditions. 

Consider now a generic density prohle, that is. 


Pi.'r) = Y. ( 38 ) 

There are physical reasons to avoid the pi term from the above expansion, namely that if 
this term is nonzero, there will be a cusp in the density at the center of the cloud. However, 
for the sake of generality, we shall not assume this term to be necessarily zero. Corresponding 
expression for F{r), and the expression for /(r) are then of the following form. 


F(r) = ^ f{r) = Y 

n=0 n=2 


(39) 


To analyse the central singularity r = 0, we now evaluate the previously dehned quantities 
? 7 (r) , j3{r) and p(r) for these density and velocity prohles, which are given by 


. 2 (2-^2 Fi ^ ^ _3r-^l SF 1 F 2 ^ SFs^ ^ ^r„l4 


V{r) = 3 + —+r[ 

-TO -TO -Co 


2 ] +-^[-..3 


.2 +-^]+0[rY- 

-CQ -Co 


(40) 
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(41) 


/9(r) = 2 + 


h ^ /2 




^ + |)] + 0[,f 


and 


A 

Fo 


Fo'f2 


+ 


Fq f2 /2-^0 / 2 -P 0 


2/21/1 

/si^i 

Af’o 

hFy /3F2 
/2ho^ / 2 -fo 

_1l + 

P 3 ^ 
-Co 


^1 

Fo‘ 


2 F 1 F 2 F 3 4 

—]+0[r] 

-TO -TO 


(42) 


respectively. Now consider the expression for 0(r) 


1 OH^ 

ei'-) = + '“' 2 ' + *=>’'1 + ; 5 (ii )72 


(43) 


where the expressions for various coefficients of d' are the following, 


'ho = 0 


(44) 


1 Js _ : rA _ , (AiAZAp 

(1 + /2/Fo)V2A2 fJ Ao f 2 ’^ i-f 2 /Fo)F^ U 2 /F 0 ) 


(45) 


1 . fs" m 2/4 f? 2F2 2F2 F,^ 

(1 + f2/Foy/^^ /2F0 + Fo' + f2 f2^ Fo^^^Fo Fo' 


3/4 3/3" sin y-Zs/Fo y + /2/F0 3sm y-^/Fo 

f2 2/2^^ (-/2/W/' (/2/F0) ^^2 (-/2/Fo)V2 
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3 + /2/-F0 w /s-^O ~ /2-^l W-^I _ 3/3 

2(1 + /2/Fo)V2^^ >^Fo /2 


(46) 


1 . 5 / 3/4 2^ _ 3/3Fi^ Fi/3^ 

(1 + /2/Fo)V2^ f,Fo h^Fo Fo^ 2 Fo3 ^ 2 / 2 ^ 0 ^ 


3/3^ , F1F4 , 3/5 3F3 3/3/4 , 3F1F2 , / 3 ^ 

8/2^0^ + 2Fo2 8^0=" ^ /2 Fo /22 Fo^ ^ 


^ 5 3 g^n-y-/ 2 /Fo _ (^f^lF,f + 4 f 2 /Fo + 3 fs _ ^ , 

Fo'^ (-/2 /Fo)3/2 /2/Fo(1 + /2/Fo)3/2 ^ ^2 Fq ^ 


3gm-y-/2/Fo _ (3 + /2/F0) -F2 Fi^ F1/3 F0/4 

^ (-/2/Fo)V2 (1 + /2/Fo)V2^^ /2 F0/2 /2' /2' 


_ Ml . ..3 gm-y-/2/Fo _ 3 + /2/F0 /3F0-/2F1 

Fo 2/2^ ^4 (-/ 2 /Fo)V 2 2(1 + / 2 /Fo)V 2 ^^ 


M M _ M . M!ii . . ^^^~V-/ 2 Mo y+ /2M0 

Fo Fo^ /2 2/2'^^ M-/2 /Fo)3/2 (/2/F0) 

^ _ 3 ^ , ^ _ M , M 4 _ 

^Fo^ Fo^ + Fo 2/2+ 2/2' 2/2=*^^ 


These coefficients completely specify the form of 0, and a is uniquely determined by the 
condition that, in the limit r —0, 0 takes a nonzero finite value. The existence of naked 


22 




singularities (black holes) is determined by the existence (absence) of real positive roots to 
the equation V{X) = 0. We can write equation (20) as (using the fact that /o = 0) 

V{X) = (1 - - X = 0 (48) 

\ X a 

where the quantity Aq is specihed by equation (31), and H{X,0) can be obtained using 
equation (14) as 

H{X,0)=X + ^ (49) 

The equation V{X) = 0, for which we need to investigate the existence of real positive 
roots, then can be reduced to a quartic of the form 

{a — l)x^ + — ©oX + = 0 (50) 

where x = \fX. The roots, positive or otherwise, of this equation can be completely 
specihed in terms of the functions ©o and Aq at the initial hypersurface. To analyse the 
existence of real positive roots in the above equation, consider a general quartic of the form 


ax'^ + Abx^ + 6cx^ + Adx + e = 0 


(51) 


with the following dehnitions, Q = ac—b"^, I = ae—ibd+Sc^, J = ace+Abed—ad"^—eb"^—c^ 
and A = I^ — 27J^. The conditions for the existence of real roots are then, A < 0 implies the 
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existence of two real and two imaginary roots and A > 0 implies that all roots are imaginary 
unless Q < 0 and (a^J — 12Q‘^) < 0 (in which case all roots are real). For the quartic under 
consideration, we have 


A = (c - 5)3A„=e„= - + (a - l)e„T (52) 

The above analysis only checks for roots to be real or imaginary; positivity or otherwise 
of the roots has to be checked separately for different cases. As explained in Section III, 
for a > 3 we only have black holes since we cannot have outgoing geodesics meeting the 
singularity in the past. Thus, for a < 3 we can have real positive roots satisfying equation 
(50). Consider first the case o; < 3. From equation (31) we then have Aq = 0, and the 
expression for A simplihes to 

A = (53) 

Since both a and ©o take only real values, it follows that A < 0 for all values of ©o, 
which ensures the existence of real roots. The roots equation V{X) = 0 gets simplihed for 
the case a < 3 as below. 


X 


3 


©0 


i.e. X = 


r Qq 12/3 


(54) 


Clearly, for the range 1 < a < 3, we always have real positive root(s) to the above 
equation, ensuring the nakedness of the singularity. However, the naked singularities in this 
case are not strong in the sense described earlier, since Aq = 0. 
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The criterion for the existence of strong curvature singularity is that a > 3, and a = 3 
is the most interesting case since a > 3 give rise to black holes. To analyze this case, using 
(50), the roots equation in this case is given as 


2a;^ + - Qqx + 


FnQo — 0 


(55) 


comparing various coefficients with the expression for a general quartic (51), we have. 


a = 2, b 




(56) 


and we have. 


A = -^0o'(40o' - 1O4 Fo3/20o + Fq^) (57) 


If the expression in parenthesis on the right side is greater than zero, we necessarily have 
two real roots to the above equation. To analyse the expression in parenthesis, consider 


g{Qo, Fo) = 400^ - 104^0^/200 + Fo^" (58) 

This equation admits roots 0i(= ^^^^[26 —5v^]) and 02(= ■^^^[26 + 5v/^]), and both 
are positive. We can easily see that the function takes negative values for the range 0i < 
©0 < ©2, and is positive outside this range. Hence the range ©o < ©i and ©o > ©2 admits 
atleast two real roots. If A > 0, we can have all the roots real if Q < 0 and the quantity 
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a^I — 12Q'^ < 0. The first condition is always satisfied, since Q = —Fo/16 (where Fq > 0), 


and the second condition is satished if 3v^(30o — < 0 , i.e. ©o < ^ 3 , say). 

Since ©3 < ©i < ©2 we cannot have real roots for the range ©i < ©o < 02 - To analyse the 
positivity of given real roots consider equation (55) which can be written as 


©0 = 



2x + Fq 

X - ^/F^ 


(59) 


The Fig. 1 explicitly shows the positivity of roots for equation (55). 



Figure 1: The graph showing existence of real positive roots of equation (55). There are no roots 
of the form x = \/To since this implies Fq = 0 which is not true. At all other points on the 
hypersurface x admits positive values with non-zero ©o (x 7 ^ 0 ). 
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It is easy to analyse from equation (15) that the regions, in the above figure, where ©o 


admits negative values necessarily give rise to shell crossings near singularity and the region 
admitting positive values of ©o is free from shell crossings for the given initial data. For a 
further discussion on shell-crossing singularities, and considerations involving more general 


class of functions than considered here, we refer to ^1 . 


To see more explicitly as to how the existence of real positive roots is related to the data 
specihed at the initial surface, (depending on the choice of F(r) and f{r)), we now analyse 
some important cases. The parameter a will decide the strength of the final singularity 
and we will be mainly concerned here for the case a = 3, i.e. strong curvature singularities 
(naked or covered). The value of a can be uniquely hxed by the limiting value of ©(r), 
which is specihed in terms of various coefficients in the expansion of T(r). In the expansion 
of T(r), if the hrst two terms Ti and T 2 vanish (Tq = 0 ), and T 3 or a higher term is the 
hrst nonvanishing term, then necessarily it is a strong curvature singularity. 


The two important specihc questions which need to be addressed in this connection are, 
(i) Given an arbitrary density distribution on the initial surface, is it possible to keep the hrst 
two terms in the expansion of T(r) zero, and later terms non-zero, by means of a suitable 
choice of a velocity prohle for the cloud, so that the resulting singularity is strong curvature 
type, and (ii) How ‘generic’ is the occurrence of the strong curvature naked singularities. 

To investigate the answer to the hrst question, we start by considering a specihc example. 
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Consider a density profile of the form 


p(r) = Po+ Pnf”, P2^0 


(60) 


n=2 


where po is the central density of the cloud. We will assume the density to be decreasing 
outwards from the center, which is a reasonable condition to take for realistic density prohles, 
and which is also a necessary condition to avoid any possible shell-crossing singularities. The 
choice of such a density distribution has some physical signihcance. Firstly, as commented 
earlier, the linear term in the density prohle is chosen to be zero, because a non-vanishing 
linear term would mean a cusp in the density at the center with a non-zero pressure gradient 
at the center. This may not be desirable on physical grounds. Next, for physically realistic 
density conhgurations, the p 2 term is generally believed to be non-zero, with p 2 < 0, for 
reasons related to the stability of the stars (see e.g. p2|, ^). The subsequent terms in 
the expansion may be vanishing or otherwise, which will not affect our considerations here. 


because as we shall see, it is the hrst non-vanishing term in the expansion which determines 
the nature of the singularity. Another useful point to note is that for a marginally bound 
collapse (/ = 0), such a density prohle necessarily results into a weak naked singularity p. 
Hence, it is of importance to know whether the collapse of the same density prohle can result 
into a strong curvature naked singularity when the collapse is not marginally bound. 


For the density prohle given above, the corresponding expression for F{r) is given by 


F(r) = For^ + Y. 




n=2 


( 61 ) 






As for the other free function /(r), the appropriate differentiability for the metric requires 


it to be at least a function. Therefore, a general choice for /(r) is, 


n=OQ 


fir) = Y. fnr 

n=2 


(62) 


We now £x /(r), which specifies the velocity profile for the cloud, by imposing the 
following two requirements: (i) the second term in the expansion of /(r), i.e. /s is zero and, 
(ii) the coefficients F 2 and /4 as given by 


F 2 




(63) 


be related as 


F 2 = u 


■{i+f2/Foy/2 


|G(- 


-/ 2 /F 0 )] 


■(i+/2/Vo)V2 -Gi-f2/Fo)] 


(64) 


Under this situation, we evalnate the expressions for r]{r) and /5(r), which tnrn ont to be 


? 7 (r) = 3 + r^[2F2] + + O^] 

ho 


(65) 


(3{r) = 2 + r‘^[2fi] + r^[^] + 0[r^] 

h 


( 66 ) 


respectively. Then, considering the general expression for 0(r), 


0 (r) 


^(r) 

^3(a-l)/2 


(67) 
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the various coefficients are evaluated to be as follows 


'ho = 0 


(68) 


'Pi ^ 0 


(69) 


= 0 


(70) 


'1'3 = 3('^)| ' 


/ 2''(1 + / 2 / Fo )‘-'= 2 


- rG(-/2/Fo)l 


(71) 


Since 'hs is the hrst nonvanishing term in the expansion of 'h(r), we have a = 3, which 
makes it to be a strong curvature singularity. Also, whenever a = 3 we know that Aq is a 
nonzero quantity, so from the roots equation V{X) = 0 we have. 




(72) 


where H{X, 0) is given by 


H{X,0)=X + 


Of 


XV 2 


(73) 


Writing again x = y/X, this becomes a quartic equation. 


+ \[Fqx^ - Qqx + Qq\[Fq = 0 


(74) 
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Depending on the choice of the values of the initial variables, that is /s, /2 and Fq, we can 


have the above quartic with or without real positive roots (Fig. 1). Hence we can have both 
the possibilities as we would desire, namely the black holes as well as naked singularities, 
depending on the choice of initial data. As we have already pointed out earlier, in either 
case this is a strong curvature singularity. 

We have illustrated above how strong curvature singularities arise as a result of collapse 
when the coefficient p 2 is the first nonvanishing term in the expansion of the density profile. 
We chose in the above the pi term in the expansion of density to be identically zero, since 
it would represent a cusp in the density at the center of the star and may be objectionable 
on physical grounds. However, it is not conclusively established that such density cusps at 
the center are ruled out by astrophysical considerations (see e.g. [^) Hence, if we assume 
that the pi term is non-vanishing and arises as a perturbation (which can possibly happen 
in the turbulent interiors of stars), in that case as well we can show that strong curvature 
singularities can arise as the final fate of collapse. For that purpose, using similar techniques 
as above, one again chooses the velocity profile suitably, imposing required constraints on 
/(r). Next, when p^ is the first non-vanishing term in the density expansion, then it is already 
known that we have a strong curvature singularity for the marginally bound collapse. This 
case is discussed earlier in detail [p, showing the formation of naked singularities and black 
holes. In a subsequent paper, the causal structure of spacetime has been analysed, using the 
dynamics of trapped surfaces, for this particular example (see Fig. 2) 
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Figure 2 


Figure 2: A plot of the apparent horizon curves (12) for marginally bound case and , p{r) = 
Po ~ ) obtained by setting po = 1 and varying rc (boundary of the cloud). For ^ < — 2 

i{= Fa/Fb^^^) the center is the hrst point to get trapped, whereas for ^ > —2 some surface in the 
interior of the star is the first one to get trapped, and the trapped region moves both inwards and 
outwards. But for all other points on the singularity curve (r > 0), ts{r) > tahii"), independent of 
the initial data. 

Interesting situation arises when all the first three terms pi,p 2 and ps in the expansion 
of density profile are zero, and the first non-vanishing term in the expansion of density is 
Pn {n > 3). We know that in the marginally bound case this situation always corresponds 
to a black hole |p. Now, considering the general non-marginally bound case, the density 
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profile is of the form 


p{r) = po+Y. Pnr"" (75) 

n=4 

The corresponding expression for F{r) can be written as 

OO 

F{r) = For^ + ^ ( 76 ) 

n=4 

Consider then the following choice of the energy function, or equivalently the particle 
velocity prohle, as given by 

n=QO 

/(r) = ^ (77) 

n=5 

We can see that for this choice of /(r) (i.e. /a = 0, fi = 0) we will have Tq, Ti and 'I /2 
identically zero, and ^3 is nonzero with 

a = 3, T3 = —(78) 

The expression for quartic remains the same equation (55), except that the functional 
form of ©0 changes now. Hence, we can have again both the possibilities, namely the 
black holes and strong curvature naked singularities depending on the choice of initial free 
functions. 
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After starting the gravitational collapse from a regular initial data, it is desirable to 
ensure that no shell-crossing singularity forms in the cloud before formation of central shell- 
focusing singularity. The necessary and sufficient condition to avoid shell-crossings in the 
cloud is given by equation (11), which can again be written as (equation (34)), 

g'(r) < 0 , s'(r) > 0 

The condition q'{r) < 0 implies that the density should decrease or remain constant as we 
move away from the center of the cloud. This condition is satished as long as the sum of all 
non-vanishing derivatives is negative (or zero), keeping the overall density positive (to satisfy 
the weak energy condition); or if all the coefficients in the density prohle are negative. We 
note that the condition for existence of strong curvature singularity derived here depends on 
the first non-vanishing term of density prohle, which is taken negative in our case, and the 
remaining terms would not contribute to the hnal outcome. Hence, rest of the terms in the 
density prohle can always be chosen to satisfy this condition by the initial data F{r) chosen 
in the analysis. Now consider the second constraint i.e., s'{r) > 0. As discussed earlier, as 
we go out from the center of the cloud, this condition implies that the contribution of kinetic 
energy to the total energy decreases. This reduces to the following condition on the initial 
data, 

F{r) < C[-rf{r)] (79) 
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Since F{r) and —rf{r) admit only positive values, C here is a non-zero positive constant. 


A lower bound on the values of C easily follows using the boundary condition at the surface 
of the star, i.e. 


C > 


M 

-rjc 


( 80 ) 


where M is the total mass of the cloud and the subscript c denotes values measured at 
the boundary of the cloud. In equation (80), M/tc and fc denote respectively the potential 
energy and the total energy of the shell at the boundary of the cloud respectively. In 
general, from equations (3) and (79) we can see that we must have C > 1 always, equality 
being satished for shells which are rest. Also, from equation (80), the equality implies that 
boundary of the star is at rest. 

Let us consider a specific example to understand the possibility of existence of regular 
initial data {F{r), f{r)) for which there are no shell-crossing singularities in the cloud till 
the hrst shell-focusing singularity forms at the center of the cloud, i.e. 


^s(O) ^ ^sc(^)|r>0 

At the center of the cloud we have r = 0,f5(0) = tsc(O). Consider a density prohle of 
the form (60). Corresponding expression for mass function is given by (61). For energy 
function /(r) of the form (62) subject to suitable constraints (/s = 0) and (64), we have 
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shown earlier that it gives rise to strong curvature naked singularity. The condition (79) for 


avoiding shell-crossing singularities reduces to 


Fo + F2r^ + - < C[-f2 - fy - fy + ■■] 

where C > 1. The above condition will hold in general, everywhere in a cloud of arbitrary 
radius, at least for that class of functions, f{r), where each coefficient in the expansion of /(r) 
satishes the above inequality with corresponding coefficients in the expansion of F{r). This 
can be easily done in this case, since /s can be given any hnite value (keeping —1 < /(r) < 0 ) 
which makes Ts 7 ^ 0 in (71), and since both /s and /2 are less then zero therefore Ts > 0, 
which is consistent with ©q > 0. Since C > 1, the argument of G{y) {y = —/ 2 /T 0 ) in (64) 
also satishes the condition 0 < y < 1. 

Now consider another example where the density prohle is given by (75). Here we make a 
choice of /(r) such that the coefficients /a and are identically zero. Clearly, the condition 
(79) is trivially satished as the corresponding coefficients in the expansion of F{r) are also 
zero. So we are free to make a suitable choice of I/ 5 I, and I/ 2 I (non-zero) for any given 
F{r) such that the required inequality (80) is satished. Hence the condition for avoiding 
shell-crossing singularities can always be satished in general without ahecting the general 
conclusions regarding the generic existence of strong curvature singularities. 

This completes our answer to the hrst question, which is in the affirmative. That is, for 
any given density distribution at the initial epoch of time, we can always keep the coefficient 
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'^3 as the first non-vanishing term in the expansion of tl/(r), by means of a suitable initial 
choice of the velocity prohle /(r) for the cloud. In other words, given any generic density 
prohle, we can always choose the rest of the initial data, that is the particle velocity prohle, 
in such a manner that the hnal fate of the collapse results in a strong curvature singularity 
which is either hidden inside a black hole, or naked, communicating with faraway observers 
in the spacetime. 

Coming to question (ii), the following consideration may provide some insight into the 
issues regarding genericity. Consider the functional form of the initial data, as we have 
discussed above. The functions F{r) and /(r) both are specihed in terms of inhnitely many 
free functions in the form of the coefficients of the expansion terms such as Fq, Fi, F 2 ,and 
/ 2 , /s,... etc. Hence, we have here a 2 ® cxd dimensional initial data space. Now keeping F{r) 
completely free and hxing hnite number of coefficients in /(r), as we have done above, to 
demonstrate the existence of strong naked singularities still leaves us with a 2®cx) dimensional 
parameter space, and hence the hnal result in terms of either a naked singularity or black hole 
is stable under a range of perturbations in the newly constrained parameter space. We would 
like to contrast this situation with the possible scenario when we could possibly have only a 
hnite dimensional initial data space which generates a strong naked singularity, as a subset of 
an inhnite dimensional initial data space. Then, since the strong curvature singularity occurs 
only for a hnite number of coefficients in the space of hnitely many coefficients characterizing 
the complete space of initial data, we could possibly say that any generic point in the 
complete space will not lie on such a hypersurface of strong singularity, and hence the data 
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generating strong singularities must be a set of “measure zero” in some sense. Such an 


argument gives an idea on the genericity of the occurrence of strong curvature singularities 
in gravitational collapse which could be either naked or covered. 


VI. Discussion and Conclusion 


We have used here collisionless fluid, that is dust with p = 0 equation of state, as the 
model to analyse the issue of the dual fate of gravitational collapse of a massive cloud, if the 
collapse started from a “physically reasonable” initial data in terms of the initial density and 
velocity prohles of the cloud. We have imposed only rather general and physically reasonable 
conditions such as the matter satisfies the weak energy condition, and that the collapsing 
dust avoids shell-crossing singularities or caustics (this is essentially because our main focus 
of interest is the nature of the shell-focusing singularity occurring at the center). For the 
sake of clarity, and to generate maximum possible physical insight, we have conhned here the 
attention to the class of Taylor expandable density and velocity prohles only and considered 
mainly the occurrence of strong curvature singularities which are naked or hidden inside 
black holes. 

An important limitation of the analysis presented here could be thought of as the choice 
of the dust equation of state. The question as to whether the introduction of pressure 
will signihcantly modify the conclusions given here needs to be looked into carefully. It 
may be pointed out, however, that one of the motivations for not considering the general 
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equation of state presently is the possibility indicated by recent work that the specihc 


properties of matter helds may not turn out to be the key factor in deciding the hnal fate 
of gravitational collapse. The indications by the work such as above are that the nature 
of the central singularity should essentially depend on the choice of initial data, and also 
there is some kind of a pattern in the behaviour of the central point (r = 0), and the other 
points on the singularity curve, regardless of the exact form of the matter used 11 , 0 , 11 . 
Again, while deciding on the issue of how signihcant (or insignihcant) the assumption of 
dust is, the arguments such as those by Hagedorn need to be considered that a very 
soft equation of state is a good approximation under extreme physical conditions of the 


advanced state of collapse. It is possible that even if we include pressures, the collisions also 
will give a contribution to the energy-momentum tensor, and will accelerate the formation 
of singularities 0 0- Hence the collisionless assumption could represent fairly general 
features of collapse and deserves serious consideration from the point of view of possible 
physical applications. From such a perspective, if we assume the formation of singularity in 
collapse, as implied by the singularity theorems as well as the physical considerations on the 
hnal fate of collapsing massive stars which have exhausted their nuclear fuel, then we have 
tried to argue here that “actual and real” strong curvature singularities (naked or covered) 
do appear in a rather “generic” way as the end product of gravitational collapse. 
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